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Abstract. This econophysics work studies the long-range Ising model of a finite 
system with TV spins and the exchange interaction and the external field H as a 
model for homogeneous credit portfolio of assets with default probability Pd and default 
correlation p&. Based on the discussion on the (J, if) phase diagram, we develop 
a perturbative calculation method for the model and obtain explicit expressions for 
Pd, pd and the normalization factor Z in terms of the model parameters N and J, H. 
The effect of the default correlation pd on the probabilities P(Nd,Pd) for Nd defaults 
and on the cumulative distribution function D(i,pd) are discussed. The latter means 
the average loss rate of the "tranche" (layered structure ) of the securities (e.g. CDO), 
which are synthesized from a pool of many assets. We show that the expected loss rate 
of the subordinated tranche decreases with pd and that of the senior tranche increases 
linearly, which are important in their pricing and ratings. 
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1. Introduction 

The statistical properties of the models for credit risks have been widely discussed in 
the past ten years from the standpoint of financial engineering |2] and econophysics 
[SI E]. In the context of econophysics, the mechanism of systemic failure in banking 
has been studied Ej- Power law behavior of the distributions of avalanches and 
several scaling laws in the context of percolation theory were found. On the other hand, 
in financial engineering, the evaluation of the effect of the correlation p between the 
rates of return of assets or between the default of assets is a hot topic and is widely 
discussed from theoretical and empirical viewpoints. Empirically, historically realized 
values of correlations and their implied values, which are estimated based on the market 
value of credit derivatives, are compared and their discrepancies, called correlation risk 
premium, attract investors' interests from the viewpoint of portfolio management [Zj. 
Theoretically, many statistical models are proposed for modeling credit risk of the pool 
of many assets |2l IHl El HDl [TTJ [T21 H31 H31 HH1 HE! - There are two categories in the models. 
The models in the first category use two state discrete variables which describe whether 
the asset is defaulted or not In the financial literature a two- valued variable 

Xi(i — 1, • • • , N) takes values and 1 depending on whether the z-th asset described by 
Xi is not defaulted or defaulted. The default probability P d is defined by the average 
number of defaulted assets Nd per an asset as Pd = <J ^ > . Here < > means the 
expectation value. Ising Spin variable S = ±1 is also used and it is related to x as 
S — 1 — 2x. Moody's Binomial (Expansion) approach |Hj|, Moody's correlated Binomial 
model |10j . Long-range Ising model ^T] are in this category. The default correlation is 
defined by the simultaneous default probability. If we denote the probability distribution 
of two asset i, j as Pij(xi,Xj), the default correlation pd is defined by 

Here, Pi(xi) and Pj(%j) are the probability distributions of Xi and Xj and they are 
calculated from the joint probability distribution Pij(xi,Xj). 

In the second category, the models adopt a continuous variable for the earning rate 
of an asset and correlation between the earning rates is introduced [HI El EH E]- On the 
assumption that the earing rates obey multivariate normal distribution with correlation 
p a , the probability Pij for the simultaneous default of the i-th and j-th assets is given 
by 

1 r*-HPi) r*~HPj) . u 2 -2p a uv + v 2 



Pi,j = , = / / exp( — )dudv. (2) 

2nJl-p 2 a J -™ J-™ 2(1 -pi) 



Here Pi and Pj are the default probabilities of the z'-th and j-th assets and $ _1 (x) is 
the inverse function of the normal distribution function. The variables u, v mean the 
earing rates of the two assets. If the random variable u (or v ) becomes lower than 
$ _1 (Pi) (resp. $ _1 (Pj)), the z-th (j-th) asset is judged to be defaulted. The correlation 
parameter p a is named as "asset correlation" and p a and pd are related via the equation 
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(HJ). The conditionally independent model [12], such as the Merton based model, the 
credit metrix model [I3J and the copula model ^3] are in the second category. 

The reason why default or asset correlations are widely discussed recently is that the 
pricing of Asset backed Security (ABS), like CDO, needs detailed information about the 
probabilities P(Nd) for defaults. Here CDO is an abbreviation for Collateralised debt 
obligation, which is a financial innovation to securitise portfolios of defaultable assets. 
The portfolio of the underlying debts (assets) collateralizes the securitites (obligations), 
CDO is a kind of ABS. Securitization by CDO, we mean to synthesize securities based 
on a pool of many assets, like loans (CLO), commercial bonds (CBO) etc. In the 
process, layered structure is introduced and securities with high priority (reliability), 
which is called senior tranche, and those with low priority (called subordinated tranche 
or equity) are synthesized. Between the senior tranche and the equity, the mezzanine 
tranche with middle priority is also synthesized. The difference between them is that 
if some of the assets in the pool are defaulted, the security with lower priority loses its 
value at first. If the rate of defaulted assets exceeds some threshold value P c , e.g. 
P c = 7% for the mezzanie tranche and P c = 10% for the senior tranche, those with 
higher priority begin to lose their values. The equity play the role of "shock absorber" . 
By the "tranche" structure, the risk of the senior tranche is reduced and investors feel 
safe about the investment. On the other hand, the interest rates of the securities of the 
equity and the mezzanine tranche are set to be higher than those of the senior tranche 
and the subordinated tranches are high-risk-high-return products. 

The default correlation becomes important when one try to estimate the expected 
loss in each tranche, which is essential in the evaluation of its price (premium). For 
example, we assume a pool of N homogeneous assets with default probability If 
there is no correlation between the defaults, P{Nd) is the binomial distribution B(N, Pd) 
and has a peak at Nj = P^N '. The standard deviation is yj 'NP^l — Pd) — y/NP^ for 
small Pd- If the threshold value P c is large enough, the upper tranche does not suffer 
from the defaults in the pool. On the other hand, in the extreme case where the default 
correlation is pa — 1, all assets behave in the same way and there are only two cases. One 
case is that all assets are not defaulted and the probability for the case is 1 — P d . The 
other case is that all assets are defaulted simultaneously and the probability is Pd- In the 
strong correlation limit (pd = 1), when there occur defaults, all assets become defaulted 
simultaneously. Both senior and subordinated tranches lose their values completely. If 
there occurs no default, both tranche does not suffer from any damage. The essential 
problem is to know the pd dependence of the probabilities P(Nd). It is important to 
estimate the expected loss rate of each tranche based on P(Nd). In addition, we should 
also study which probabilistic model is good or useful in order to describe the behaviors 
of the assets. 

This paper deals with these problems. The organization of the paper is as follows. 
In section |21 we study the (J, H) phase diagram of finite size long-range Ising model 
and show that the assets begin to be correlated in the "Two Peak" Phase in the (J, H) 
plane. The realistic magnitude of the default correlation ranges from 1 % to several % 
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[2], only the Two Peak Phase is interesting from the financial engineering viewpoint. 
Section El is devoted to the calculation of the important parameter Pd and pd in terms of 
J, H and N. Here, we develop a perturbation method which is based on the discussions 
in section |21 Up to zero-th order in the perturbation theory, P(N d ) is expressed as the 
superposition of two binomial distributions, corresponding to the two peaks of P(Nd). 
The developed method and obtained relations are useful when one apply the long-range 
Ising model to the evaluation and hedging of the securities with tranches. In section 0J 
we study the pd dependence of P(Na, Pd) and of the expected loss rates of the tranches. 
For the latter purpose, we introduce the cumulative distribution D(i,p d ) and discuss 
that they are directly related with the average loss rates of tranches. As the correlation 
becomes strong (with fixed default probability Pd), the left peak becomes taller and 
moves towards to the origin (Nd = 0). The right peak also becomes taller and shifts 
to Nd = N. Its area approaches to Pd as pd comes close to 1. These behaviors are 
different from those of the binomial expansion approach, where P{Nd) has only one 
peak and its shape becomes broader as pd increases. We then discuss the pd dependence 
of D(i,pd). D(i) for large i increases linearly with pd and the senior tranche cannot 
avoid the default damage of the assets pool, even when we set P c to be large. This 
crucial behavior of the long-range Ising model has been pointed out previously fT] , we 
have clarified the importance in the evaluation of the tranches. Section |5] is dedicated 
to concluding remarks and future problems. We discuss the usefulness of the long-range 
Ising model from the viewpoint of financial engineering. 

2. Model and Phase Diagram in (J, H) plane 

We use Ising Spin variables S\, 5*2, ■ • • , Sn = ±1 which represent states of assets in the 
reference pool. Here Si — — 1 indicates default of i-th asset and Si = 1 means that 
the i-th asset is not defaulted. We denote the number of S — ±1 spins by N±, so 
the number of defaulted assets Nd is iV_ . The probability distribution for the states of 
the assets is assumed to be described by the following canonical distribution with the 
long-range Ising model of a finite system with N spins and the exchange interaction 
■jf and the external filed H, which are measured in units of Boltzmann constant times 
temperature. 

1 ( J N \ 

P(S 1 ,S 2 ,-..,S N )= exp — £ SiS j + Hj2Si). (3) 

We do not omit the % — j terms in the Hamiltonian for later convenience. As is well- 
known, the exchange interaction —j^SiSj controls the strength of the correlation between 
Si and Sj and the external field H favors one of the two spin states. In the actual case 
where the spin variable represents the states of the assets, the default probability Pd is 
at most a few percent and almost all assets are not defaulted (S — 1). The sign of the 
external field H is set to be if > 0. 

The reason to choose the long-range Ising model is that it gives the default 
distribution P{Nd) directly. In ^T], another motivation for the long-range Ising model 
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has been discussed and their conclusion is that the model is the most natural choice from 
the viewpoint of the Maximum Statistical Entropy principle. The two parameters J and 
H are introduced as Lagrange multipliers which ensure that the default probability and 
the default correlation of the model are pd and pd- From the economical viewpoint, 
we can interpret the model as a kind of factor model. Here, the term 'factor' means 
the systematic risk factor or the state of the business cycle In a boom, we have 
fewer defaults than in a recession. We denote the state of the business cycle as H' 
and assume that the defaults of the assets are independent from each other, conditional 
on the realization of the systematic factor H' . The joint probabilities for the assets 
Si, S2, ■ ■ ■ , Sn and the business cycle variable H' is assumed to be written as 

1 f N 



PfactASi, S 2 , ■ ■ ■ , S N , H') = —— exp \H' Si x P(H'). (4) 

^n{J) V i= i ) 

Here, the random variable H' obeys the probability density function P{H') and the 
denominator Z^(J) is the normalization term. Condition on the realization H' = H, 
the each asset state becomes independent from each other and the default probability 
Pd is given as 

-2H 

P d = Prob(S, = -1) = ^. (5) 

The default probability Pd is a decreasing function of H and H for a boom (recession) 
is large (small). In order to derive the long-range Ising model starting from the above 
factor mode, we assume that H' obeys the standard normal distribution with mean H 
and variance J/N. 

P(H') = -L= exp(- (iJ/ ~5 )2 ) (6) 
1 ; V2^J 1 2J/N ' 1 ; 

By averaging over the possible realization of H' weighted with the above P(H'), we 
obtain the expression for the long-range Ising model. 

/oo roo i / N \ 

PfacUS!, S 2 ,---, S N , H')dH' = / —— exp [H"£Si)x P{H')dH' 
-00 J— 00 ZjAr(J) \ ,1 / 



1 



N 



ZMBj^J^W + HX*). it) 

The validity of the Maximum Statistical Entropy principle or the factor model with the 
normally distributed business factor H' should be checked by the comparison with other 
more reliable models. 

The Hamiltonian of the long-range Ising model depends on the spin variables only 
through the combination of the magnetization M = J2iL% Si- There is a simple relation 
between iV_ = Nj and M as M = N + — iV_ = N — 2Nd, the default number distribution 
function P{N d ) is 

exp(^N 2 + HN) f2 J _ 

P(N d ) = z MH) L NC Nd exp N 2 d - (2 J + 2H)N 
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Figure 1. Phase diagram in (J, H) plane. For large H and small J, P(Nd) has a 
single peak at Nd — PdN. We call the region as "Single Peak" Phase. For small H 
and large J, there are two peaks in P(Nd) and we call the region "Two Peak" Phase. 

The phase boundary is depicted with the broken line (- -). The solid line ( ) 

corresponds to a constant Pd line. It starts at (J, H) = (0, He), where pd — 0. In the 
J — > oo limit, the line approaches H = ^k- asymptotically and pd — ► 1. 



The default probability Pd is defined by the expectation value of N a as 

(») 

Here < > is the expectation value with the probability distribution (|SJ). For J = 0, 
the probability distribution (jHJ) becomes that of the binomial distribution B(N, Pd) and 
there is a relation between H and as 

H = ~\ l0g {l^P d ) f ° r J = °- (10) 
We denote this value of H as He- On the other hand, for J — > +oo limit, there are 
only 2 configurations with all spins up or all spins down that have nonzero probabilities. 
The probabilities are 

1 e~ 2HN 

P(M,---,i)= 1 + e - 2HN P(-i,-i,---,-i) = 1 + e _ 2g „ - (ii) 

From the relation P(— 1, —1, • • • , —1) = Pd, one obtains the following relation between 
H and Pd for J = oo as 
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For general J, H and N, it is difficult to obtain Pj. However, for large enough N, by 
changing variable from iV_ to n_ = in eq.fJHJ), we can estimate Pd =< n_ > by the 
saddle point approximation. The saddle point equation is 

AJn--(2.J+2H) 

n - = e 4Jn--(2J+2H) _j_ l ( 12 ) 

Of course, by changing variable from n_ to the magnetization per spin m — 1 — 2n_, 
the saddle point equation is transformed into the famous self-consistent equation of the 
magnetization m = tanh( Jm + H) ^7]. Depending on the values of the parameters 
J,H(> 0), there are two cases. For large H and small J, the equation ()12)) has only 
one solution n_*. We call this region in the (J, H) plain as "One Peak" Phase, because 
the probability distribution P(Nd) has a single peak at iVn_*. Pd is almost the same 
with in the One Peak Phase. For small H and large J, the equation f)12j) has three 
solutions, two are at maxima n._i* < n_2* and one is at minimum. We call the region in 
the (J, H) plane as "Two Peak" Phase, as the reader may easily anticipate the reason. 
In the case, there is no simple relation between Pd and the solutions n_i*, n_*2- If H > 
is large, the solution n_i* is almost the same with Pd- However, when the correlation 
pd is large, the strength of H is of the order of O(^) and we cannot neglect the second 
peak n_* 2 - In the case, n^u < Pd and the average value of n_i* and n_ 2 * with P(ra_*i) 
and P{n^2) corresponds to the value of Pd- For example, when pd = 1 and J = oo, 
the average value of n_i* = and n_ 2 * = 1 with probabilities eq.(JTTJ) is equal to Pa- 
in figure ^ we summarize the situation. The solid curve (— — ) in the (J, H) plane 
corresponds to the constant Pd line. The dotted line (- -) is the "phase transition" 
line between the One-Peak Phase and the Two Peak Phase. In the remainder of the 
section, we study the correlation pd in the (J, H) phase diagram. We will see that pd is 
almost zero in the One Peak Phase. Only in the Two Peak Phase pd can take nonzero 
value. 

We discuss the default correlation pd and recall the definition (JTJ). In order to obtain 
P(Si, S2), in equation (jH}, we take the trace over S3, S4, ■ ■ ■ , Sn- 

N 

P(S 1 ,S 2 ) = H £ P(S 1 ,S 2 ,---,S N ). (13) 

j= 3s j= ±l 

The trace over S3, ■ ■ ■ , Sn is replaced by the summation over iV_ of iV" = N — 2 Spins. 
We obtain 

p(s " S2) = OT) exp (^ S2 + i/(Sl + S!) ) x 

N" 

J2 A r,,C A r_e ( ^ (5l+52)(Ar ''~ 2Ar - )+ ^ (Ar ''~ 2Ar - )2+H(Ar ''" 2Ar - )) . (14) 
7V_=0 

If the system size N is large, the summation over iV_ is replaced by the estimation at 
the saddle points. In the One-Peak region, the saddle point is at = Pd and P(S\, S 2 ) 
is given by 

P(Sx, S 2 ) oc exp(-^5i5 2 + H(St + S 2 )) exp(-J(5' 1 + S 2 ) x 2n_*)- (15) 
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Figure 2. Second peak contribution in pd- Along the line with ea. (|17fl . we plot 

the approximated estimation for pd- The lower solid line ( ) shows the data from 

ea. (fT5|) and the data from ea. lfTHjl are depicted with upper dotted curve ( ). 

In the Two Peak region, two saddle points contribute to the summation and P{Si, S 2 ) 
is estimated as 

P(S 1 , S 2 ) oc exp{^3 1 3 2 + H(S 1 + S 2 )) x 

[ e -J(Si+S a )x2n_i. x p(n_ u ) + e -J(S^S 2 )x2n^, x p( n _^).] ( 16 ) 

Here P(n_i*) and P(n_2*) are the probabilities for the two peaks and n__ 2 *- In 
the One Peak phase, the constant Pd line in the (J, H) plane is almost given by the 
following relation between J and H 

H = tanh -1 (l - 2P d ) - J(l - 2P d ). (17) 

We calculate the default correlation p d with equations (fT3j) and (fTfij) on the above 
approximate constant Pd line. About the two saddle points n_i*,n_2* and their 
probabilities P(n_i*), P(n— 2*), we take them the values at J = 00 and pd = 1. We set 
n_i* = 0, n„ 2 * = 1 and P(n_i*) = 1 — Pd, P(n^ 2 *) = Pd- We set Pd = 0.01 and we plot 
Pd vs J in figure El We see that the correlation with equation (fTKJ) . which is plotted with 
solid line (— — ), does not become large even in the Two Peak Region. On the other 

hand, p d with equation (flfi|) . which is depicted with dotted line ( ) becomes large 

in the Two Peak Region. We see that the existence of the second peak in P(A^) plays 
a crucial role in the emergence of correlation in the long-range Ising model. 

3. Perturbative Calculation and Second Peak Contribution 

In this section, we try to calculate several quantities of interest of the long-range Ising 
model. In particular, we obtain the expressions for P d and p d in terms of the model 
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parameter J, H and N. In addition, we also obtain the expression for the probability 
(or weight) of the second peak P a u, which means that almost all assets are defaulted 
11 . The probability P a u plays a crucial role when one discuss the evaluation of the 
tranche. 

When one calculate Pd, Pd, one way is to calculate < Si > and < SiSj >. Here, 
we calculate the moment of with the probability distribution eq.(jHJ). The default 
probability is then given by Pd = < ^ > . About the default correlation p d , we start from 
the following relation. 

P(-l,-l)-P(-l) 2 < S1S2 >-<S 1 xS 2 > 
Pd P(-l)(l-P(-l)) (l-<5i>)(l+<5i>) 1 j 

The magnetization M = N < S\ > and Nd is related as M = N — 2Nd and 
iV(iV — 1) < S1S2 >=< M 2 > —N, we obtain the following expression 

< + ^« Nj >-N<N d » 
Pd <N d >(N-<N d >) ■ 

In order to calculate the moment < N l d > with eq.©, the quadratic term exp(^iVj) 
prevents us from taking summation over iV^. As we have noted previously, the 
distribution with J = is binomial distribution and taking summation over Nd is 
easy. In addition, the Pd is at most a few percent and the distribution P{Nd) have a 
peak very close to N d = (and the second peak at N d — N in the Two Peak Phase). 
We expand the quadratic term as 

-p(^Hgi(^. (20, 

and perform the calculation of the moment < N l d > perturbatively. The expansion is 
about jj-Nj, which is evaluated as j^-Nj ~ 2JNPJ. In the actual risk portfolio problem, 
Pd is at most 2 ~ 3% and the system size iV is several hundred, the perturbative 
approximation is considered to be applicable. We also note that, in the Two Peak 
Phase, the above expansion should be carried out also at Nd = N. 

In order to perform the calculation in more concrete manner, we use variables N± 
and start from the following expression for the Hamiltonian. 

-H = —M 2 + HM 
2N 

J 2 7 

= —N 2 + HN - (2H + 2J)N_ + — N 2 (21) 

j 2 7 

= N 2 -HN - (2.1 - 2K)N+ + — N 2 (22) 

2N v N v ' 

In the vicinity of AL = 0, we denote P(A^_) as P_ (AL) and we can expand the quadratic 

term in eq. ()21)l . Likewise, in the vicinity of AL = N(N + = 0), we call P(A^_) as P + (A^_) 

and it can also be expanded in AL. 

P_(AL) = - N C N _e HN e^ 2H+2J ^ + ^ N - 

z 

= ± N C N _e HN e-( 2H +^ N - x £ h 2 4 N ~) k (23) 
Z ^ Q k\ N 
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p +(AL) = l CN _ N _ e -HN e -(2J-*HXN-N-)+%(N-N-)> 
Zj 

= \ N C N - N _e- HN e-^ J -^ N - N -^ x £ L&n - AL) 2 ) fc .(24) 
Z is the normalization constant to ensure that En.=o p ( n -) = !■ To the 

zero-th order 

perturbation approximation P_(AL) and P + (N_) are binomial distributions and P(AL) 
is given by the superposition of these distributions. We summarize the situation as 

( PJN) (0<AL<L) 
P(N) = { K ' K ~ ~ ' (25) 
[P+(AL) (0 < N-N- < N-L). 

Here L is set to be at the middle of the interval [0, N]. 

The moment < N l _ > is calculated with the following equation. 

N L N 

<N l _>= P(N^)N l _= P-(N^)N[+ J2 P+(N-)N l _. (26) 

7V_=0 N-=0 N-=L 

The summation over AL is from to L, however P_(AL) damps rapidly in AL, it is 
not so bad to change the range from [0, L\ to [0,A/]. About P + (AL) the range of AL 
is [L,N]. We change variable from AL to AL = N — AL and denote the probability 
distribution P+(AL — N — AL) also as P+(AL). 

P + (AL) = }- N C N+ e- HN e-^ 2J - 2H ^ + ^ N ^ . (27) 

Zj 

P + (N + ) also damps rapidly in AL, we will change the summation range from [0, N — L) 
to [0, AT]. < N l _ > is then calculated perturbatively as 

N N 

<N l _>= P-(N-)N l _ + E P +{ N +){ N ~ N +f 

JV_=0 N + =0 

= \e HN E NC N ^ H ^ N -±^Nl) k N l _ 

^ V_=0 fc=0 ^ v 

+ -e"™ }Z N C N+ e-^-^+ £ h^Nl)\N - N + )K (28) 

^ N + =0 k=0 K - iV 

The normalization constant Z is calculated as 

N oo i r> t 

Z = e HN N C N _e^ H ^ N -Y^4 N ~ )k 

Ar_=0 fc=0 fc - ^ V 

N OO i r\ J 

+ e' HN E NC N+ e~^-^Y^4 Nl)k 

N+=0 k=0 K - iV 

= Z_ + Z+. (29) 

In equation (J29j) . we denote the two terms as Z±, which come from the summation over 
AL and AL. 

In the above calculation, moments of the binomial distribution appears frequently. 



We introduce the following unnormalized binomial moments 



N l ± 



N 



N l ± 



E NC N± e^N l ± . (30) 

Af±=0 
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Figure 3. The default correlation pd on a constant Pd line. Pd = 0.01 and we plot 
Pd versus J. H is set to be the value which realize Pd = 0.01 for N = 50. The solid 

line ( ) shows the result from the perturbative calculation up to second order in 

^ and the line with + symbols shows numerically calculated exact data. 



The parameters f3± are defined as (5± = —2{J^H). Calculations of 
straightforward. The zero-th moment [1±] is given by 

[U] = (!+<*)". 



N l ± 



is 



(31) 



The /-th moment 
repeatedly. 

r jvi' 



N l ± 



is then obtained by differentiating [1±] with respect to (3± 



d l 



d(5, 



(32) 



We show the results for the first 6 moments, which are necessary for the second order 
perturbative calculation. 



[Ni 


= (} + e^] 


N 

) (N lX ±) 


(33) 


N 2 ± 


= (l+e fe " 


) N (N lX± + N 2 xl) 


(34) 


Ni 


= (l + e^ 


) N (N lX ± + 3N 2 x 2 ± + N 3 xl) 


(35) 


Ni 


= (l + e^ 


f (N lX ± + 7N 2 x 2 ± + 6N 3 x% + N 4 x 4 ± ) 


(36) 


Ni 


= (l + e^ 


f (N lX ± + 15N 2 x 2 ± + 25N 3 x 3 ± + 10^4 + N 5 x 5 ±) 


(37) 


Ni 


= (i+e & ; 


f (N lX ± + 31N 2 xl 





+ 90N 3 x 3 ± + 65N 4 xi + 15N 5 x 5 ± + N 6 x 6 ± ), (38) 
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N l ± 



is 



where iVj = r^zW an d x ± = ]T~kt ■ I* 1 general, the l-th binomial moment 
calculated as 

N l ± ] = (l + e^) N j2y l ±k N k x k ± , (39) 
fc=i 

where the coefficients y l ± k for N k x± is calculated with the following recursive relations. 

y l ±k = y l ±U + ky 1 ^ (40) 

and with the conditions y l ±k = for k > I and y\ = 1. 

With these preparations, we are ready to write down the results. The perturbative 
calculation of the normalization constant Z is given as 

z = Z^ + Z + 
The moment < N l _ > is given by 



T 2k+h 



Z f^LkWN) 

oo -r / n T\ k I 



+ ^' HN Et,(4) E iC m (-l) m N l - m [Nl k+m }. (42) 



z fe fo^Viv; m=0 

Putting these results into Pd = <7 ^ d> and eq. (|TU|) . the expressions for P d and in 
terms of model parameters N, J, H are obtained. In addition, the weight of the second 
peak P a u, that is the probability of almost all assets are defaulted, is estimated as 

Pall = (43) 

As we have noted previously, the zero-th order approximation means that we express 
P(Nd) as a superposition of two binomial distributions. In the case, the results for Z 
and Pd, Pd can be written down in the following simple expressions. 

Z° = e HN (l + e /3 -) N (l + 1 N ) (44) 



x_ + 7^(1 — X H 

1 7 ^ 



p° d = ::s +; (45) 

(l-(x_+x+)) 2 (46) 



The subscript (0) indicates the zero-th order perturbation results. In figure El we shows 
the result for p\ along the constant P d line. P d is set to be Pd = 0.01 and with solid line 
(— — ) we show the data from the above perturbative calculation up to second order in 
The line with + symbols depicts the numerical data. The two lines coincide well 
and the match is very good as long as Pd is set to be small. Figure 0] is the 3-dimensional 
plot of the data in (J, H, p d ) space, pd begins to be large in the Two Peak region and 
its rapid growth is well captured by the above perturbative calculation. 



Evaluation of Tranche in Securitization 



13 




Figure 4. 3-dimensional plot of p d in (J, H, pd) space. The solid line ( ) shows the 

result from the perturbative calculation up to second order in 2j£ and the line with + 
symbols depicts numerically calculated exact data. The conditions are the same with 
those in figure |3| 



4. Effect of p d on P(N d ) and on average loss rates of tranches 

We would like to discuss the effect of default correlation p d on the probabilities P{N d ) 
and on the tranche synthesized from the pool of the homogeneous assets. In order to 
discuss the latter case, we introduce the cumulative distribution function D(i, p d ), which 
is directly related with the average loss rate of the tranche. The default rate P d and the 
system size N is fixed. When we show numerical data, we set N = 100 and P d = 0.05. 

At first, we discuss the former case, the effect on the probabilities P(N d ,p d ). Here 
we write down their p d dependence explicitly. The p d = case is easy and P(N d , 0) is 
only the binomial distribution B(N,P d ). It has a single peak at NP d and the width 
is roughly 2^/NP d . For N d = 0, P(0, 0) = (1 — P d ) N and as p d becomes large, 
P(0,p d ) approaches 1 — P d . In figure we plot P(N d ,p d ) versus p d for N d < 10. 
The system size N = 100 and P d = 0.05. P(0,p d ) grows monotonically as p d grows. 
For 1 < N d < 5 = NP d , P(N d ,p d ) at first increases and then decreases as a function 
of p d . On the other hand, for N d > 5 = NP d , P(N d ,p d ) decreases with p d . P d is 
small and P(N d ,0) damps rapidly in N d for N d > 5, P(N d ,p d ) is almost zero for 
any p d , which holds for 10 < N d < 90. Figure |H1 depicts the plots of P(N d , p d ) for 
90 < N d < 100. P(N,0) = P£ ~ and P(N, 1) = P d , P(100,p d ) grows monotonically 
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Figure 5. Plot of P(N d , p d ) vs p d . < N d < 10 and N = 100, P d = 0.05. 



from to P d — 0.05. For 90 < N d < 99, P(N d , p d ) is upward convex with respect to 
Pd- The area of the second peak becomes greater with the increase of pd and P(N d , pd) 
increases for N d ~ N. As pd becomes large, the width of the second peak becomes 
narrow and P(N d , p d ) for N d ^ N decreases. On the other hand, P(N, p d ) increases 
monotonically to P d . 

To sum up, for p d = 0, P(N d ,0) is B(N,P d ) and it has a single peak. The 
width of the peak is order 2y/NP d and it is small for small Pd- As pd grows, the 
system is in the Two Peak Phase. At the zero-th order perturbative approximation, 
P(N-) = P-(N-)+P + (N + ) is a superposition of two binomial distributions. P_(N d ) has 
a peak at Nd < NPd and is approximately obeys B(N, On the other hand, P + (Nd) 
is B(N, 1 — x + ) and has a peak at Nd ~ N. The increase in p d accompanies the increase 
in J, however the change of H is not so large and it decreases slightly (See Figure HJ). 
The first peak position of P(N d ), which is governed by /3_ = —2(J + H), moves towards 
to Nd = as J increases. The first peak becomes narrower and higher with the left 
slide and only P(0, pd) grows monotonically. For < N d < NP d , P(N d , p d ) is upward 
convex with respect to pd- P(Nd,pd) for N d > NPd damps with pd monotonically. On 
the other hand, the second peak position, which is governed by /?+ = — 2( J — H), shifts 
towards to N d = N. The area of the second peak, which is calculated as in eq.()43|). 
approaches Pd and the width becomes narrow. P(N, pd) increases monotonically to Pd 
with p d and P(N d , p d ) near N d = N is upward convex. 

We would like to discuss the above effect on the tranche of securities synthesized 
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Figure 6. Plot of P(N d , p d ) vs p d . 90 < N d < 100 and N = 100, P d = 0.05. 



from the homogeneous assets pool with parameters P d: p d - For the purpose, it is useful 
to introduce the cumulative distribution function D{i,p), which is defined as 

N 

D(i,p d )= ]T P(N d ,p d ). (48) 

N d =i 

From the definition D(0,p d ) = 1 and D(l,p d ) = 1 — P(0,p d ) is the probability of the 
occurrence of default. We explain the relation between D(i,p d ) and the evaluation of 
the tranche briefly. 

The tranche for the interval implies that if the number of default N d is below 
% (N d < i), the tranche does not suffer from any damage. However, if N d exceeds or 
becomes equal to i (N d > i), it begins to lose its value. The value of the tranche is 
A = j — i + 1 in units of the number of assets (we assume that the values of all assets in 
the pool are equal.) and if defaults with i < N d < j occurs, it loses (N d — i + units. 
When N d exceeds j(N d > j), the tranche lose its value completely. The expected loss 
rate of the tranche is calculated as 

l ( 3 N \ 

E(i\j) = -r(j:P(k,p d )(k-(i-l)) + A £ P(k,p d )\. (49) 
^ \k=i k=j+i ) 

The first terms comes from the partial damage in the tranche {% < N d < j) and the 
second term implies the contribution from its complete loss of the tranche (N d > j + 1). 
E(i\j) are directly related with the price of the tranche (premium), which can be 
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observed in the market. For j = i, we denote E(i\i) as E(i) and call it as the expected 
loss rate at the z-th tranche. It is related with the cumulative distribution D(i,p) as 

N 

E{i) = E{i\i) = E P(k, Pi ) = D{i, Pd ). (50) 

k=i 

E(i) is useful, because we can reconstruct E(i\j) as a sum of E{k) as 

k=i 

The proof of the relation is straightforward. 

j j n j j N 

E^(0= E(E^p,)) = E(En^) + E p (k,p d )) 

l=i l=i k=l l=i k=l k=j+l 

j N 

= J2 p (k,P d )(k-(i-l))+A £ P(k) = AE(i\j) (52) 

We note that, if we set i = 1 and j = N in equation (|5ip. we obtain 

1 * 

£(l|JV) = -5>(*)=-P* (53) 
fc=i 

Here we use the relation £7(1 |iV) = Pd, which is intuitively clear and can be proved as in 
the equation (J52|) . From the second equality in eq.(jHSJ) that the average of the expected 
loss rate at each tranche is P d , tranches look like to "share Pd between them" or "toss 
Pd to other tranches". 

Now we discuss the effect of the default correlation p d on E(i) = D(i, p d ). From the 
definition D(i, p) = Y,k=i P{k, Pd), we can understand the pd dependence easily from the 
previous discussions on P(N d ,Pd)- In figure [3 we show D{i, p d ) for 91 < i < 100. The 
area of the second peak increases monotonically to Pd as we increase pd, the cumulative 
distributions D{i, pd) also grow up to Pd- P(Nd,pd) is almost zero for 11 < N d < 90, 
D(i,pd) for 11 < i < 90 behaves in the same way with these D(i,pd) for i = 91. If 
i becomes small, we see the contribution from the first peak in P(Nd,pd)- For small 
i, the damps of P{N d , Pd) for N d ^ with respect to p d dominates the contribution 
from the second peak. P(0,pd) increases monotonically and D(l,p d ) = 1 — P(0,pd) 
decreases with p d . D{i,p d ) for i > 2 also decreases as in figure |Hl which shows D{i,p d ) 
for 1 < i < 10. These behaviors reflect the left shift and the width tinning of the first 
peak. For the intermediate value of i, the p d dependence of D{i,p d ) is not monotonous. 
In figure 01 we depict D(i,p d ) for 9 < % < 11. Along with the shape change of the first 
peak with p d , D(i,p d ) at first decrease. Then, the contribution from the second peak 
dominates the decrease of the first peak contribution and D(i,p d ) begins to increase. 
D(i,pd) is downward convex with respect to pd for the interval of i. 

We note that the ranges where D{i, p d ) is downward convex, D{i,p d ) decreases 
monotonically, or D(i,pd) increases monotonically depends on the parameters N,Pd- 
The above discussions may not hold for other values of N and Pd- In particular the 
range of the downward convex region, if we set Pd = 0.01, we observe that it shift to 
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Figure 7. Plot of D(i,p d ) vs p d . 91 < i < 100. 
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the left. The positions of the boundaries between the regions are important from the 
viewpoints of risk management and rating of the securities, we should note this point. 



0.1 i 1 1 1 r 




02 I I r r~~~ I I I I I I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Pd 



Figure 9. Plot of D(i, pd) vs pd- 9 < i < 11. Solid line i = 9, broken line i = 10 and 
dotted line i = 11. 

From these observations, we summarize the pd dependence of E(i) = D{i,p d ). 

• Senior tranche, the range with i,j are set to be large, suffers from the 
default correlation seriously. E(i) for the range increase linearly with pd- It is 
approximately given by 

E(i) = Pd xP d (54) 

If pd change from to 1, E(i) change from to Pd and the evaluations of the 
securities decrease almost linearly. 

• Equity or subordinate tranche, the range is [1, j] and j is small. In the range 
1 < k < j, D{k,pd) damps monotonically with pd and E{l\j) also damps. The 
increase in pd causes the increase of the evaluation of the tranche. 

• Mezzanine tranche, whose range is between the equity and the senior tranche. In 
the range, the behaviors of D(k,pa) depends on Pd and the system size N. In the 
above mentioned iV = 100, Pd = 0.05 case, D(k,pd) has downward convex shape 
for some intermediate values of k. 
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5. Concluding Remarks and Future Problems 

In this paper, we have studied the long-range Ising model as a model for a pool of N 
homogeneous assets with default probability P d and default correlation p d . We have 
studied in the (J, H) plane, the behavior of Pd and p&. There are two phases in the 
(J, H) plane. In the One Peak Phase, the probabilities P{Nd) have a single peak at 
Nd ~ NPd- The correlation is almost zero in the phase. In the Two Peak Phase, there 
are two peaks in P{Nd) and pd can take large value. The first peak is closer to origin 
than NPd and its area is larger than 1 — P^. The second peak is at about N and its area 
is less than Pd- The parameters J, H should be chosen in the Two peak phase, if the 
model intends to describe the portfolio with some default correlation between the assets. 
We have developed the perturbative method and expressed P(N d ) as a superposition of 
two binomial distributions with the above two peaks at zero-th order. We have obtained 
the closed form expression for Pd, pd and the weights for the second peak P a u, which 
means the probability that almost all assets are defaulted. These expressions are in 
good agreement with numerically calculated values and give an efficient method for the 
actual application of the long-range Ising model. Otherwise, for Pd and pd, it is difficult 
to know the parameters J and H and the long range Ising model is hard-to-use as a 
model for homogeneous credit risk portfolio. 

Furthermore, we have studied the pd dependence of P(Nd, Pd) and the cumulative 
distribution D(i,p d ). P(N d ,0) is binomial distribution B(N,P d ) and it has a peak at 
NPd- As we increase pd from to 1, the profile of P(Nd,pd) changes from One peak 
shape to Two Peak shape. The first peak shifts to the left and its shape becomes higher 
and narrower. The second peak's area increases and it shifts to the right with the 
decrease of its width. At pd = 1, P{Nd, 1) has two thin peaks at N d = and N d = N 
and the probabilities are P(0, 1) = 1 — Pd and P(N, 1) = Pd- Other probabilities are 
zero. The cumulative distribution functions D(i, pd) correspond to the average loss rates 
of the z-th tranche. About the senior tranche, the range of the tranche is large. 
As pd increases, D(i,pd) increase almost linearly with pd like D(i,pd) — Pd X Pd- The 
average loss rate of the senior tranche is given as a sum of D(k,p d ) in the range 
the expectation value of the loss rate of the senior tranche also increases as pd x Pd- 
The price of the tranche is based on the average loss rate, the value of the senior tranche 
decreases with pd- The range of the equity, the subordinated tranche, is near the origin 
and the D(i,p d )s decrease monotonically. The average loss rate of the equity decreases 
with pd and the price of the equity increase with pd- The mezzanine tranche is between 
the equity and the senior tranche. The profile of D(i,pd) in the range depends on the 
model parameters pd,Pd and N. In the text example, D(i,p d ) has a downward convex 
shape in some region. If the mezzanine range is chosen to lie in the region, the 
average loss rate also behaves similarly. However, other probabilistic model for a pool 
of assets, e.g. the copula model [H], suggest upward convex shape for the average loss 
of the mezzanine tranche. The discrepancy comes from the difference of the shapes of 
P(Nd,Pd)- The more complete comparison between the probabilistic models for a pool 
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of correlated assets should be done. 

As concluding remarks, we comment on the usage of the long-range Ising model 
and related future problems. As a statistical model for an ensemble of many assets, 
the long-range Ising model is an attractive one from the viewpoint of physicists. Its 
phase diagram and phase transitions are throughly studied and its analytic calculation 
method, like Hubbard-Stratanovich transformation, guides us how to make theoretically 
tractable models. On the other hand, from the viewpoint of financial engineers, the long- 
range Ising model is not so convenient. One reason is that the model parameters J, H 
are not directly related with the observed data P d and p d (or p a ). Other statistical 
models incorporate these parameters as a model parameters. For example, the Moody's 
correlated binomial model gives P(N d ) as a function of P d and p d explicitly. When one 
uses Ising model, it is necessary to know the parameters J, H which correspond to P d , p d . 
The definition of P d and p d include the moments < N l _ > or < Si > and < SiSj >, it is 
necessary to take the trace Tr. The long-range Ising model has the advantage that the 
trace Tr is reduced to the summation over the total magnetization M = Y,f=i Sj and 
the calculation is not so heavy task. Even so, this one step spoils the usefulness of the 
model. We have obtained a closed form expressions for P d and p d and try to circumvent 
the step. The computational time to obtain J, H for given P d , p d is reduced much and 
the failing of the model are partially overcome. 

In order to apply the long-range Ising model to the evaluation of the tranche [i.j] 
in more realistic situation, the assumption of homogeneity of the assets pool should be 
weakened. One step toward the direction is to introduce multi sectors and assume the 
homogeneity only in each sector. We label each sector by / = 1, 2, • • • K and i-th sector 
contains N 1 assets. In the 7-th sector, the default rate is Pj and the default correlation 
is p d . Between different sectors, say between I-th and J-th sector, the default correlation 
is p 1 /. We use Ising Spin variables Sj to represent the states of the j-th asset in the 7-th 
sector, the generalized long-range Ising model Hamiltonian for the probabilities P(Sj) 
is 

- « = £ 4v M - + H ' M - + vwm M ' Mj (55) 

As in the homogeneous model, the Hamiltonian depends on Sj only through the 
magnetization of the J-th sector Mj = J2fli Sj. When we set j^- = ^p^p? = jf 
and H 1 = H, the model reduces to the homogeneous model with N = Ylf=i N 1 , J, H. 
The problem is to get the relation between p d , p 1 / , P d and J 1 , J IJ , H 1 , N 1 . In order 
to accomplish the task, the phase diagram in J IJ , J 1 , H 1 and the profile P(iVj) should 
be cleared and it is left for future analysis. Furthermore, for more complex situation 
where i-th asset has default probability P d and the default correlation between i-th and 
j-th asset is p d , the model Hamiltonian becomes that of the random Ising spin systems. 
The exchange interaction and the external field Hi should be connected to P l d and 
p d , which is also left for future problem. Other step is to discard the Ising model and 
adopt other probabilistic models. One possibility is the Moody's correlated binomial 
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model, which uses two state variables x = 0, 1 for the state of an asset and incorporates 
Pd and pd directly in the model parameters. Its generalization to the multi-sector case 
and more complex situations is an interesting problem. Other possibility is to introduce 
simplified version of the long-range Ising model. We use two state variable X; for the 
state of the z-th asset. The number of defaults Nj, is expressed as = J2jL\Xj. The 
probabilities P(xi, x 2 , ■ ■ ■ , xn) is given as 

P(x u x 2 ,---,x N ) = Q.-a)p N *(l-p) N - N * + a6 NiNd . (56) 

Instead of the superposition of two binomial distributions, we use a5N,N d for the second 
peak. The first peak is B(N,p) and the parameters a,p are related with the default 
probability Pd as Pd = a+(l—a)p. This probabilities P(xi, x 2 , • • • , xn) is more tractable 
than the original probability distribution Q and the generalizations to more complex 
situations may be carried out easily. 
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